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Abstract. In this paper we investigate Gaussian queues in the light-traffic and in the 
heavy-traffic regime. Let = {Qx : * — 0} denote a stationary buffer content 

process for a fluid queue fed by the centered Gaussian process X = {X(t) : t £ R} with 
stationary increments, X(0) = 0, continuous sample paths and variance function <r 2 (-). 
The system is drained with a constant rate c > 0, so that for any t > 0, 

Q x ) (t) = sup (X(t)-X(s)-c(t-s)). 

— oo <s<f 

We study = {Q^x (t) : t > 0} in the regimes c — > (heavy traffic) and c — > oo 
(light traffic). We show for both limiting regimes that, under mild regularity conditions 
on a, there exists a normalizing function 5(c) such that Qx\d(c)-)/a(3(c)) converges to 
( ) in C[0, oo), where Bh is a fractional Brownian motion with suitably chosen Hurst 
parameter H. 



1. Introduction 

A substantial research effort has been devoted to the analysis of queues with Gaussian 



input, often also called Gaussian queues [lCj, [ill, Il2j ]. The interest in this model can 



be explained from the fact that the Gaussian input model is highly flexible in terms of 
incorporating a broad set of correlation structures and, at the same time, adequately 



approximates various real-life systems. A key result in this area is [181 ]. where it is shown 
that large aggregates of Internet sources converge to a fractional Brownian motion (being 
a specific Gaussian process). 

The setting considered in this paper is that of a centered Gaussian process X = {X(t) : t G 
M} with stationary increments, X(0) = 0, continuous sample paths and variance function 
<7 2 (-), equipped with a deterministic, linear drift with rate c > 0, reflected at 0: 

Q ( x\t)= sup (X(t)-X(s)-c{t-s)). 

—oo<s<t 

The resulting stationary workload process can be regarded as a queue The objective 
of the paper is to study = {Qx CO : i > 0} in the limiting regimes c — >• (heavy 
traffic) and c — > oo (light traffic). 
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Under mild conditions on the variance function c 2 (-), Q$ is a properly denned, almost 
surely (a.s.) finite stochastic process. However, if c — > 0, then Q x (t) grows to infinity (in 
a distributional sense), for any t > 0. The branch of queueing theory investigating how 
fast Q$ grows to infinity (as c — > 0) is commonly referred to as the domain of heavy- 
traffic approximations. In many situations this regime allows manageable expressions for 
performance metrics that are, under 'normal' load conditions, highly complex or even 
intractable, see for instance the seminal paper by Kingman [9] on the classical single- 
server queue. Since then, a similar approach has been followed in various other settings, 
see, e.g., @, 0, 15, 17, 19] and many other papers. 



Analogously, one can ask what happens in the light-traffic regime, i.e., c — > oo; then 

(c) 

evidently Q x decreases to zero. So far, hardly any attention has been paid to the light- 
traffic and heavy-traffic regimes for Gaussian queues. An exception is Debicki and Mandjes 
0], where the focus is on a special family of Gaussian processes, in a specific heavy-traffic 

(c) 

setting. The primary contribution of the present paper concerns the analysis of Q x under 
both limiting regimes, for quite a broad class of Gaussian input processes X. 

We now give a somewhat more detailed introduction to the material presented in this 
paper. It is well known that under the assumption that cr(-) varies regularly at infinity 
with parameter a € (0, 1), for any function 5 such that 5(c) — > oo as c — > 0, there is 
convergence to fractional Brownian motion in the heavy-traffic regime: 

(1) ^§^4^(.),asc^0. 

a(5(c)) 

We shall show that an analogous statement holds in the light-traffic regime, that is, if a(-) 
varies regularly at zero with parameter A G (0,1) (i.e., x h-> o~(1/x) varies regularly at 
infinity with parameter —A), then for any function 5 such that 5(c) — > as c — > oo, 

( 2 ) ^^4s A (-),a8c^oo. 

a(6(c)) 

Assuming that X satisfies some minor additional conditions, both ([1]) and ([2]) apply in 
C(R), the space of all continuous functions on R. 

Our paper shows that the statements ([T]) and ([2]), which relate to the input processes, carry 

(c) 

over to the corresponding stationary buffer content processes Q y x - That is, we identify, 
under specific conditions, a function 5(-) such that 

Qg(-),asc^0 



and 



Qx (<?(c)0 4^(1), 
a(5(c)) 

Q [ x } (Ml < n" > 



a(5(c)) 



Qb'A'), as oo, 



both in the space C[0,oo) of all continuous functions on [0,oo). 

This paper is organized as follows. In lSection 21 we introduce the notation and give some 
preliminaries. ISubsection 3.11 presents the results for the heavy-traffic regime, whereas 
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ISubsection 3.21 covers the light-traffic regime. We give the proofs of the main theorems, 
(i.e.. IThcorcm ll and IThcorcm 2j) in lSection 41 



2. Preliminaries 

In this paper we use the following notation. By id : R — > R we shall denote the identity 
operator on R, that is, id(t) = t for every t € R. We write f(x) ~ g{x) as x — > xq G [0, oo] 
when lim x _j, Xo f(x)/g(x) = 1. Let MY^ip) and g&'f'^X) denote the class of regularly 
varying functions at infinity with parameter a and at zero with parameter A, respectively. 
That is, for a non-negative measurable functions /, g on [0, 00), / G (ffi'¥ Q)0 (ot) if for all 
t > 0, f(tx)/f(x) -)• t a as x -)• 00; 5 G (A) if for all i > 0, g{tx)/g{x) ->• i A as x ->■ 0. 

2.1. Spaces of continuous functions. We refer to Billingsley [3J for the details of this 
subsection. For any T > 0, let C[— T, T] be the space of all continuous functions / : 
[— T, T] — > R. Equip C[— T, T] with the topology of uniform convergence, i.e., the topology 
generated by the norm ||/||[_t,t] := su Pte[-T,T] 1/(01 under which C[—T,T] is a separable 
Banach space. Therefore, by Prokhorov's theorem, weak convergence of random elements 
{X®} of C[-T,T] as c 00 is implied by convergence of finite-dimensional distributions 
and tightness. A family {X( c )} in C[—T,T] is tight if and only if for each positive s, there 
exists an a and cq such that 

(3) P(jX (c) (0)| > aj < e, for all c> c ; 

and, for any rj > 0, 

/ 



(4) lim lim sup I 

C~ >"0 c— >oo 



sup 

. |t-s|<C 
\s,t€[-T,T] 



x ( c \t) - X^ (s) 



> n\ =0. 



For notational convenience, we leave out the requirement s,t G [— T, T] explicitly in the 
remainder of this paper. 

Finally, let C(R) be the space of all functions / : R — > R such that /ir-j^l G C[— T, T] 
for all T > 0. The above definitions extend in an obvious way to C[0, T], C[0, 00) and 
convergence as c — » 0. 

For 7 > 0, let il 7 be the space of all continuous functions / : R — > R such that 
lim^ioo f(t)/(l + \t\ J ) = 0. Equip f2 7 with the topology generated by the norm ||/||m := 
sup tGK |/(t)|/(l + |f| 7 ) under which 7 is a separable Banach space, so that Prokhorov's 
theorem applies. The following property can be found in 0, Lemma 3] or [J, Lemma 4]. 

Proposition 1. Let a family of random elements {X^} on f2 7 be given. Suppose that the 
image of {X^} under the projection mapping pt '■ ^ — > C[—T,T] is tight in C[—T,T] 
for all T > 0. Then {X^} is tight in f2 7 if and only if for any r\ > 0, 

/ |X (c) (t)| \ 
(5) lim lim sup P sup j— — > 77 =0. 

T-s>oo c _>oo \|t|>T 1 + It / 
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2.2. Fluid Queues. Let Q$ = {Q^x(i) : t > 0} denote a stationary buffer content 
process for a fluid queue fed by a centered Gaussian process X = {X(t) : f £ 1} with 
stationary increments, X(0) = 0, continuous sample paths and variance function 0" 2 (-). 
The system is drained with a constant rate c > 0, so that for any t > 0, 

Q { x\t)= sup (X(t)-X(s)-c(t-s)). 

—oo<s<t 

Additionally, an equivalent representation for Q^{t) holds 16, p. 375]: 

(6) Q ( x\t) = Q ( x\0) + X(t) -ct + max (0, sup f-Q^(0) - (X(s) - cs))) . 

V 0<s<t v ' J 

Throughout the paper we say that X satisfies: 

C: if er 2 (t)| log |i|| 1+e has a finite limit as t — > 0, for some e > 0; 

RV : if a G W (A), for A G (0, 1); 

RVoo: if er € MV^a), for a G (0, 1); 

HT: if both C and RVoo are satisfied. 

LT: if both RVo and RVoo are satisfied. 

Remark 1. In our setting (X has stationary increments), the assumption that X is continu- 
ous is equivalent to the convergence of Dudley integral, see lSubscction 231 This is immedi- 
ately implied by condition C, see [l|, Thm. 1.4]. However, the real importance of condition 
C lies in the fact that if in addition X satisfies RVoo, then X also belongs to il 7 , for every 
7 > a. This is pointed out in ISubsection 3.l1 Finally, note that C is met under RVo. 
Indeed, since a G @y (\), then t \-> a(l/t) belongs to ^^(-A), thus cr 2 (l/t)t A -> 
as t — > oo. Equivalently, a 2 (t)t~ x — > as t — > 0, implying lim t ^of 2 (*)| log |£|| 1+E = 0, for 
any fixed e > 0. Furthermore, RVqo implies that X{t)/t — > a.s., for t — >• =too, so that 
is a properly defined stochastic process for any c > 0, see [T], Lemma 3]. Lastly, the 
assumption that X has continuous sample paths implies that a is continuous. 

Due to the stationarity of increments, all finite-dimensional distributions of X are specified 
by the variance function, since we have 

(7) Cov(X(i),X( S )) = \ (a\s) + a 2 (t) - a\\t - s\)) . 

Recall that by Bh = {Bait) : t G M} we denote fractional Brownian motion with Hurst 
parameter H G (0, 1), that is, a centered Gaussian process with stationary increments, 
continuous sample paths, Bh(0) = and covariance function 

(8) Cov(B H (t),B H (s)) = 1 {\s\ 2H + \t\ 2H -\t- s\ 2H ) . 

As mentioned in the introduction, if c — > 0, then, for any t, Q^(t) — > oo a.s., which is 
called the heavy-traffic regime. On the other hand, if c — > oo, then Q^(t) — > a.s., which 
is called the light-traffic regime. 
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2.3. Metric entropy. For any Tel define the semimetric 



d(t, s) := y/E\X(t) - X(s)\ 2 = a(\t - s\), t, s G T. 

We say that S C T is a $-net in T with respect to the semimetric d, if for any t € T 
there exists an s £ S such that d(t, s) < ~d. The metric entropy Hd(T, i?) is defined as 
logNd(T, •&), where N^(T, •&) denotes the minimal number of points in a i?-net in T with 
respect to d. Later on we use the following proposition, see 0, Thm. 1.3.3] and [2, Corollary 
1.3.4], respectively. 

Proposition 2. There exists a universal constant K such that for a d-compact set T 

rdiam(T)/2 



/ \ /•diam(T)/2 

E ( supX(t) I < K I v/e d (T,<!?)d??. 
\tef J Jo 



and for all £ > 
/ 

E 



sup \X(t) - X(s)\ | < K y/Ud (T, #) di?. 



. (s,t)eTxT 
V d(s,t)<C 



The quantity J °° (T, ??) di? is called the Dudley integral. 

3. Main Results 

In this section we formulate the result for the heavy-traffic and light-traffic regime, respec- 
tively. It is emphasized that these results are highly symmetric. Let us first introduce a 
function 5, such that for every c > 

By the continuity of a, we can choose 5 as 5(c) = inf{x > : xjo(x) = 1/c}. From the 
definition of 5 it follows that 5 G MVq(1/(ol - 1)) under RVoo and 5 £ Mf oc (l/(X - 1)) 
under RVq. 



3.1. Heavy-traffic Regime. In the heavy-traffic regime we are interested in the analysis 

(c) 

of Q x as c — > 0, under the ass 
follows from 0, Thms. 5 and 6]. 



(c) 

of Q x as c — > 0, under the assumption that X satisfies HT. The following statement 



Proposition 3. If X satisfies HT, then 
o-(<5(c)) 

in C(M) and O 7 , /or any 7 > a. 

In fact, |Proposition~3 holds for any function 5(c) such that 5(c) — > 00 as c — ?■ 0. Condition 
C (which is one of the requirements of HT) plays a crucial role in proving tightness both 
in C[—T,T], for some T > 0, and in $7 7 . 

Combining [Proposition 3 with the definition of 5 leads to the following statement. 
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Corollary 1. If X satisfies HT, then 
X(5(c)-) -ctf(c) id(-) d 



a(S(c)) 

in 



B a (-) -id(-) as 0, 



Now we are in the position to present the main result of this subsection. 
Theorem 1. If X satisfies HT, then 



(10) "X™' 4QW(0 asc^O, 



tr(ff(c)) 
m C[0, oo). 

We postpone the proof of IThcorcm ll to lScction~4l 



Remark 2. Theorem [T] extends the findings of [8|, Theorem 3.2] where, under the heavy- 
traffic regime, the weak convergence in C[0,oo) of Q { $ (5(c)-) /a(5(c)) as c -> was obtained 
for the class of input processes having differentiable sample paths a.s., i.e., of the form 
X(t) = f Z(s)ds, where {Z(s) : s > 0} is a stationary centered Gaussian process whose 
variance function satisfies specific regularity conditions. 

(c) 

3.2. Light-traffic Regime. In the light-traffic regime we analyze the convergence of Q x 
as c — > oo, under the assumption that X satisfies LT. We begin by stating the counterpart 



of Proposition 3 



Proposition 4. If X satisfies RVo, then 

X(5(c)-) d 

777 ^ B M")> as c ^ oo, 
a(5(c)) 

in C(R). //, moreover, X satisfies LT, t/ien i/ie convergence also holds in ft 7 , for any 
7 > max{A, a}. 



Analogously to Proposition 3 Proposition 4 holds for any function 5(c) such that 5(c) — > 
as c — >• oo. As in the heavy-traffic case, combining Proposition 4| with the definition of 5 



leads to the counterpart of Corollary 1 



Corollary 2. If X satisfies RVo, then 
X(5(c)-)-c5(c)id(-) d 



a(5(c)) 

in C(R). 



B\(-) — id(-) as c-> oo, 



The main result of this subsection is now stated as follows. 
Theorem 2. If X satisfies LT, then 
nu Qx\Hc» d n (i) n 

in C[0, oo). 
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We postpone the proof of Proposition 4 and ITheorem 21 to ISection 4l 



Remark 3. The assumption LT excludes the class of input processes of the structure 
X(t) = Jq Z(s)ds, with {Z(s) : s > 0} being a centered stationary Gaussian process with 
continuous sample paths a.s. (since A = 1 in this case). In [8|, Theorem 4.1] it was shown 
that, for this class of Gaussian processes, Q£ (0)/er(<5(c)) does not converge weakly to 
Q^(0) as c ->■ oo. 

4. Proofs 

In this section we prove our results, but we start by presenting an auxiliary result. 

Lemma 1. If X satisfies LT, then for any e > 0, there exist constants C, a > 0, such that 
for all x < a and t > 0, 

cr(x) ~ 1 t u t > 1, 
where I := min{A — e, a + e} and u := max{a + e, A + e}. 

Proof Take any e > 0, then because a £ i^^o (-*•)> there exists an a < 1 such that 

(12) < 2t x ~ e , for all x < a and ix < a. 
a(x) 

Moreover, there exists a constant K\ such that cr(x) > K\x x+e for all x < a. 
Because a 6 ^"^^(a), there exist constants A,K2 > such that <r(x) < K2X a+e for all 
x > A. Because a is continuous, we can in fact find a K% such that cr(x) < -K2# a+e for all 
x > a. Therefore 

cr(tx) i^2(ix) a+e , n ^ „ 

-V^ < A rr =: Kt a+e x a ~ x , for all x < a and tx > a. 
cr(x) ~ KiX A+e 

Note that, if q — A > 0, then we have 

(13) < Ka a+e t a+e , for all x < a and tx > a. 
a(x) 

If a — A < 0, then 

(14) < Ka a ~ x t x+e , for all x < a and tx > a. 
a(x) 

Combining ()12j) (|14)) . we conclude that there exists a constant C > 0, such that 

< Cmax|t A ~ e ,t Q+e ,t A+e ), for all x < a and all t > 0. 
<r(x) I J 

□ 

In what follows, we will use the following notation. Let 



xW(t);= M 

<t(<5(c)) 
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and denote the variance of by (a^) 2 , that is, 
a[5[c)) 



Proof of Proposition 4| We begin by showing the convergence in C(R). To this end, we 
need to show the convergence in C[—T,T] for any fixed T > 0. 

Convergence in C[—T,T]: From the fact that a G ^^o(X), it is immediate that the finite- 
dimensional distributions of converge in distribution to B\ as c — > oo, cf. ©-([8]), 
which also implies d3|). Therefore, the weak convergence of X^ in C[-T,T] follows after 
showing Q. 

By the Uniform Convergence Theorem, see 0, Thm. 1.5.2], for any t G (0, £], we have 
cr( c )(i) < 2( x . Thus, Proposition 2 yields, for some universal constant K > 0, 



sup 

>-*l<C 





> r? j < IP 


sup 


X^(t) - X {c \s) 


> 






\cr(c) (|s-t|)<2C A 







I< c )(t)-l( c '(s) 



< - E f sup 

^ \(j( c )(|s-t|)<2C A 

<_/ JM')([-T,T\,0)d&, 
V Jo v 



where H( c )([— T, T], •) is the metric entropy induced by a^ c \ 

By Potter's bound |4|, Thm. 1.5.6] for any e, £ > 0, e < A and t G (0, £] and sufficiently 
large c (corresponding to small 5(c)), we have cr( c )(t) < 2t x ~ e . Hence 



M( c )([-r,T],^)<M d -([-r,T],0 



where d is a semimetric such that d(s, t) = \t — s\ x e . The inverse of x \— > x x e is given by 
x i — y x l ^ x ~ e \ so that 

m-TW) < + 1) < ciog (i 



for some constant C > and $ > small. It follows that 

^HM([-T,T],tf) d#<VC~£ C ^log (j)dti = 2y/cj* 
Summarizing, we have 

2KVC f°° Vlog? 



limsupP sup 

c->oo \|s-i|<C 



we obtain ((31) by letting £ — >• 0. 



#2 



\/logi? 

I? 2 



di9; 
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Convergence in f2 7 : To show the convergence in $7 7 , we need to verify ([5|). Observe that 

l^Esu Pte[eJ>e3+1] |X( c )(t)| 



sup 



\xvm > \ < if. 



1 ^ E\X^(ei)\ 2 ^ Esu Pte[eJ - eJ+1] X^(t) 



< 



77 ^ 1 + eJT 77 ^ 1 + e-^ 

■h(k) + h{k). 



I\{k) and /2(&) are dealt with separately. According to lLemma 1\ for large c (that is, small 
8(c)), we have 



a( c \t) <Cx 



r t<i, 
t u t > l, 



where £ and it can be chosen such that £,u < j. Therefore, 

v ; ~ 77 ^ 1 + eii ~ 77 ^ 1 + 
' j=* ' i=* 

and the resulting upper bound tends to zero as k — > oo. 

Now focus on l2(k). For some universal constant K > and because of the stationarity of 
the increments of X, [Proposition 2] yields that !%$) is majorized by 



,,, x oo y V /H( c )([ei,ei+ 1 ],7?) dtf ^ oo J y%( c )([0,e>'(e-l)],#) 



77 ^ 1 + e-JT 77 ^ 1 + e-JT 

j=k j=k 

We will estimate the integrals under the sum by splitting the integration area into i9 < 1 
and t? > 1. 

Observe that, for some constants C\, C% > (that is, not depending on j), 



£ y/m(°)([0,ej(e-l)],#)d# < £ Jlog ( ^jf + l) dT? 



./Ci+j 



/■oo 

/ v^e- M d79 = C 2 e^. 

JO 

Recall that £ < so that 

^/H( c )([0,eJ(e-l)],7?)dtf 



oo / yM" U u,e^e-i;j,i/;ui/ oo ,y 

lim V ^ : < lim V 2 . = 0. 

fc^oo ^ 1 + e J7 ~ 77 fc^oo ^ 1 + &7T 
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So it remains to show the analogous statement for the integration interval [l,oo). Using a 
similar argumentation as the one above, one can show that 

^/e( c )([0,ei(e-l)],tf) d& < C 3 e uj , 

for some constant C3 > 0, from which the claim is readily obtained. □ 

Since the proof of lTheorem H is analogous to the proof of lTheorem 21 we choose to focus 
on the light-traffic case only. 

Proof of lTheorem 21 The proof consists of three steps: convergence of the one-dimensional 
distributions, the finite-dimensional distributions, and a tightness argument. 

Step 1: Convergence of one- dimensional distributions. In this step we show that, for a 
fixed t > 0, 

^ 4 Qg> (*),»<: -00. 

(c) 

Since Q x is stationary, it is enough to show the above convergence for t = only. Observe 
that, due to the time-reversibility property of Gaussian processes, 

<#(0) = sup (X(t) - ct) = sup (X(5(c)t) - c8(c)t) . 

t>0 t>0 



Upon combining Corollary 2 with the continuous mapping theorem, for each T > 0, 



/ X(S(c)t) - c5(c)t\ d 
SU P 777 ^^ ^ SU P ( S a(*) ~ t), as c -> 00. 

te[o,T] V °"l d ( c )j / te[o,T] 

Thus it suffices to show that 

(15) lim lim sup P ( sup ( ^(c)t) -^c)t \ \ 

T^oo c ^oo \t>T V °"(^( c )) / J 

for any 77 > 0. Recall the definition of X^ c \ so that 

/ ( X(5{c)t) - c5{c)t \ \ ( \X^{t)\ \ 

P SUp -— ■ > 7] < P SUp — — > 1 , 

\t>T\ cr{S{c)) J J y t > T n + t J 
where we used Q. Proposition 4 implies that the family {X( c )} is tight in , for some 



7 < 1. Now (I15p follows from Proposition 1 



Step 2: Convergence of finite- dimensional distributions. The argumentation of this step is 
analogous to Step 1. First note that for any t{ > 0, rji > and Sj < ti, where i = 1, . . . ,n, 
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for any n E N, it follows that 

F { a(5(c)) >Vi,i = h...,n) 

= pf sup ^Wc^-xw-cc^-.n > \ 

<rf sup (W.)-^)-*fe-»)\ >% , = 1| ^\ 

+ v p (sup > , 

Now the same procedure can be followed as in Step 1. 

S'tep 5: Tightness in C[0,T]. In this step, for any T > 0, we show the tightness of 
{Q^ (5(c)-) /a (5(c))} in C[0,T]. Given that we have established Step 2 already, ([3]) holds 
so we are left with proving Q, with s,t E [0,T]; the remainder of the proof is devoted to 
settling this claim. 

Stationarity of implies that {Q^ (5(c)t) — Q^(5(c)s) : t > s} is distributed as 

{Q^(5(c)(t-s))-Q^(0):t>s}, 
so that it suffices to prove (jlj) for s = only. Furthermore, cf. ([6]), 
Q^(5(c)t) - Q5?(0) < 2 sup |X(<5(c)t) - c5(c)t\ . 



sup 

o<t<c 



o<t<c 



From Corollary 2 it follows that 



\X(5(c)t) - c5(c)t\ d 
sup sup \B\(t) — t\ , as c — >• oo. 

o<t<c °"(<H C )) o<t<c 



Now notice that for £ < r//4, by the self-similarity of S 



A; 



sup \B\(t) - t\ > \ ) < 2P ( sup B\(t) > ^C X ) 

3<t<C 2 / \0<t<l 4 / 



\o<t<C 

Now it is straightforward to conclude that the last expression tends to zero as £ — > 0. □ 
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